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Introduction 

The  use  of  sulfur  mustard  [SM,  bis-(2-chloroethyl)  sulfide] 
as  a  chemical  warfare  agent  dates  back  to  World  War  I. 
Unfortunately,  there  is  no  current  treatment,  and  relatively  little 
is  known  about  which  proteins  play  a  major  role  in  SM-induced 
toxicity.  Indeed,  a  milieu  of  signaling  molecules  respond  to  SM, 
triggering  many  pathways  including  apoptosis,  inflammation, 
and  DNA  damage  (7).  Recent  reports  have  shown  that  SM 
exposure  activates  specific  proteins,  such  as  p53,  NF-/cB,  Bcl- 
2,  and  others  (7).  While  these  studies  are  effective  at  probing 
individual  pathways,  they  do  not  put  into  context  the  global 
changes  that  are  occurring  in  response  to  SM  and  how  these 
many  signaling  pathways  intersect;  thus,  they  provide  only  a 
limited  snapshot  of  SM-induced  toxicity. 

Over  the  past  few  years,  innovations  in  systems  biology  and 
biotechnology  have  led  to  important  advances  in  our  under¬ 
standing  of  biological  phenomena.  This  is  due  in  part  to  better 
instrumentation  and  software  development.  In  particular,  the 
field  of  proteomics  has  been  accelerated  by  more  accurate  and 
high-throughput  mass  spectrometry  instrumentation.  Protein 
samples  that  once  took  hours  to  evaluate  can  now  be  analyzed 
in  a  matter  of  seconds. 

Technologies  that  allow  for  rapid,  sensitive,  and  specific 
analysis  of  the  cellular  proteome  have  already  proven  effective 
at  answering  a  range  of  biological  questions,  from  targeted 
studies  to  global  analyses.  For  example,  metabolic  labeling 
approaches  such  as  stable  isotope  labeling  with  amino  acids  in 
cell  culture  (SILAC)  (2)  have  proven  valuable  for  studying  cell 
culture  models  of  disease  and  cellular  perturbations.  SILAC  has 
been  used  to  study  tumor  metastasis  (3,  4),  focal  adhesion- 
associated  proteins,  growth  factor  signaling,  and  insulin  regula¬ 
tion  (5).  More  specifically,  SILAC  has  been  employed  to  study 
post-translational  modifications  (PTMs),  including  methylation 
(5)  and  phosphorylation  (6). 

The  SILAC  method  employs  12C14N-  and  13C15N-labeled 
amino  acids  added  directly  to  the  culture  medium,  resulting  in 
isotopically  “light”  and  “heavy”  cell  populations,  respectively. 
Protein  samples  collected  from  control  (light-labeled)  and 
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experimental  (heavy-labeled)  cells  can  then  be  mixed  in  equal 
ratios,  digested  with  trypsin,  and  analyzed  by  high-resolution 
mass  spectrometry.  The  corresponding  light  and  heavy  peptides 
from  the  same  protein  will  coelute  during  chromatography  into 
the  mass  spectrometer.  Because  these  peptide  ions  have  the  same 
ionization  efficiency,  relative  quantitative  information  is  obtained 
for  each  protein  by  comparing  the  signal  intensities  of  the  light 
and  heavy  peptides  in  the  precursor  ion  scan  (MS  scan). 
Sequence  information  is  then  gathered  via  peptide  fragmentation 
in  the  fragment  ion  scan  (MS/MS  scan).  Thus,  the  incorporation 
of  stable  isotopes  allows  for  a  quantitative  comparison  between 
the  two  samples.  Furthermore,  SILAC  can  be  coupled  to 
additional  PTM  enrichment  methods  to  provide  a  more  special¬ 
ized  understanding  of  PTM-regulated  cellular  processes. 

Among  these  PTM  enrichment  methods  is  immobilized  metal 
affinity  chromatography  (IMAC)  (7)  for  the  isolation  of  phos- 
phorylated  peptides.  In  short,  IMAC  employs  positively  charged 
metal  ions  that  will  bind  to  negatively  charged  phosphopeptides 
from  a  digested  lysate.  After  an  incubation  period  to  allow 
binding,  samples  are  washed  to  reduce  nonspecificity,  and  then, 
phosphopeptides  are  eluted  for  subsequent  analysis.  Recently, 
this  approach  has  been  combined  with  mass  spectrometry  to 
provide  a  relatively  quick  and  effective  means  of  enriching  and 
identifying  thousands  of  phosphopeptides  from  biological 
samples  (6). 

To  date,  neither  SILAC  nor  IMAC  has  been  applied  to  the 
study  of  SM,  and  the  use  of  other  proteomics  methods  relating 
to  SM  has  been  limited.  In  a  recent  proteomics  study  (8)  of 
cells  treated  with  radiolabeled  SM,  two-dimensional  gel  analysis 
demonstrated  that  SM  forms  adducts  with  cytoskeletal  proteins. 
While  this  work  provided  valuable  information  on  specific 
targets  of  SM,  it  offered  a  relatively  small  amount  of  data  given 
all  of  the  downstream  signaling  that  occurs  in  response  to  SM 
exposure.  Alternatively,  SILAC  has  the  ability  to  identify  and 
quantify  hundreds  of  proteins  in  SM-treated  cells,  leading  to 
unique  insight  into  the  biological  response  to  SM-induced 
toxicity  as  well  as  offering  a  significant  increase  in  the  number 
of  potential  therapeutic  targets. 

Because  a  large  amount  of  cellular  signaling  activity  is 
mediated  by  PTMs  such  as  phosphorylation,  here,  we  present  a 
large-scale  quantitative  analysis  of  the  phosphoproteome  using 
SILAC  and  IMAC  in  a  human  keratinocyte  model  of  SM- 
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induced  toxicity  [immortalized  human  keratinocyte  line  (HaCaT) 
cells  (9,  70)].  Furthermore,  we  show  that  the  levels  of  numerous 
protein  phosphorylation  sites  are  altered  in  response  to  SM  and 
that  many  of  these  resulting  changes  in  SM-exposed  cells  were 
previously  unknown. 


Experimental  Procedures 

Cell  Culture,  Metabolic  Labeling,  SM  Exposure,  and  Lysate 
Preparation.  For  SILAC  labeling,  HaCaT  cells  were  cultured  in 
arginine  (Arg)-  and  lysine  (Lys)-deficient  Dulbecco’s  modified 
Eagle’s  medium  (DMEM)  supplemented  with  10%  dialyzed  fetal 
bovine  serum  (FBS),  sodium  pyruvate,  glutamine,  and  gentamycin/ 
amphoteracin  and  maintained  in  15  cm  dishes  at  37  °C  for  at  least 
five  population  doublings.  Media  were  also  supplemented  with 
either  12C614N4-Arg  and  12C614N2-Lys  (“light”  condition)  or  13C615N4- 
Arg  and  13C61;,N2-Lys  (“heavy”  condition).  Arg  was  added  to 
DMEM  at  one-third  of  the  manufacturer’s  recommended  amount. 
Once  cells  were  at  approximately  80%  confluency,  light-conditioned 
cells  were  treated  with  vehicle  (control),  and  heavy-conditioned 
cells  were  exposed  to  200  /<M  SM  (experimental),  a  concentration 
that  induces  inflammatory  cell  signaling  (10,  11)  for  20  min  in  a 
certified  chemical  surety  fume  hood.  For  whole-cell  lysate  protein 
preparation,  plates  were  washed  2x  in  ice-cold  PBS  and  collected 
in  lysis  buffer  containing  20  mM  Tris-HCl  (pH  7.5),  1%  Triton 
X-100,  1  mM  EDTA,  150  mM  NaCl,  and  protease/phosphatase 
inhibitor  tablets  (Roche,  Indianapolis,  IN).  Cells  were  homogenized 
with  glass  beads  and  a  mechanical  shaker  for  30  s  and  ultracen- 
trifuged  at  20000g  for  10  min  at  4  °C.  Pellets  were  discarded,  and 
supernatant  protein  concentrations  were  determined  by  BCA  assay 
(Pierce,  Rockford,  IL). 

SDS-PAGE  Fractionation,  In-Gel  Digestion,  and  IMAC 
Enrichment  of  Phosphopeptides.  Two  milligrams  of  total  protein 
(1  mg  of  light  lysate  mixed  with  1  mg  of  heavy  lysate)  was  resolved 
by  SDS-PAGE  and  briefly  stained  with  colloidal  Coomassie  blue. 
After  it  was  destained,  the  gel  lane  was  excised  into  six  regions, 
and  each  region  was  cut  into  1  mm  cubes.  Processing  for  in-gel 
Uypsin  digestion  was  performed  as  described  (72).  Following  trypsin 
digestion,  IMAC  (7)  enrichment  of  phosphopeptides  was  carried 
out  as  described  (13). 

LC-MS/MS  Analysis.  LC-MS/MS  experiments  were  performed 
on  a  LTQ  Orbitrap  hybrid  mass  spectrometer  (Thermo  Electron, 
San  Jose,  CA)  as  described  (14). 

Phosphopeptide  Identification  and  Quantification.  MS/MS 
spectra  were  searched  via  the  SEQUEST  algorithm  (75)  against  a 
composite  database  containing  the  human  International  Protein 
Index  (IPI)  protein  database  and  its  reversed  complement  (16). 
Parameters  allowed  for  up  to  three  missed  cleavages,  a  mass 
tolerance  of  ±100  ppm,  and  a  static  modification  of  57.0215  Da 
for  alkylation  on  cysteine.  Dynamic  modifications  included  79.9663 
Da  for  phosphorylation  on  serine,  threonine,  and  tyrosine;  15.9949 
Da  for  oxidation  on  methionine;  and  10.0083  and  8.0142  Da  on 
Arg  and  Lys,  respectively,  to  account  for  the  use  of  isotopically 
encoded  amino  acids.  Search  results  were  filtered  to  include  <1% 
matches  to  reverse  sequences  by  restricting  the  mass  tolerance 
window  and  setting  thresholds  for  Xcorr  and  dCn'  (16).  We  used 
the  Ascore  algorithm  for  phosphorylation  site  localization  (77). 
Unique  phosphorylation  sites  were  determined  by  utilizing  the  best 
Ascore  for  each  phosphorylation  site  to  calculate  the  minimum 
number  of  possible  phosphorylated  residues  on  each  individual 
protein.  Phosphorylation  site  quantitation  was  performed  using  Vista 
(18).  In  cases  where  the  same  phosphorylation  site  was  quantified 
multiple  times,  the  best  example  in  terms  of  Vista  confidence  score 
and  minimum  S/N  was  used  (14).  This  qualitative  and  quantitative 
workflow,  including  SILAC  labeling,  IMAC  enrichment,  MS 
analysis,  and  peptide  identification/quantitation,  was  carried  out  in 
triplicate  (n  =  3  biological  replicates). 

Western  Blotting.  HaCaT  cells  were  cultured  as  described  above. 
A  time— course  study  was  conducted  with  cells  either  untreated  or 
SM-exposed  (200  fiM)  for  0,  10,  20,  30,  and  60  min.  Lysates  were 


_ Control _  SM-treated 

Time  (min)  0  10  20  30  60  0  10  20  30  60 

—  — •  — -  Phospho-p38 
Total  p38 

Figure  1.  Confirmation  of  p53  and  p38  phosphorylation  following  SM 
exposure.  Western  blot  analysis  of  HaCaT  cells  that  were  unexposed 
(control,  left  panel)  or  exposed  (SM-treated,  right  panel)  to  SM  for  0, 
10,  20,  30,  and  60  min.  Results  showed  low  levels  of  phosphorylation 
for  both  p53  and  p38  in  the  control  samples  but  a  steady  increase  in 
phosphorylated  p53  and  relatively  intense  p38  phosphorylation  by  20 
min  postexposure.  Total  levels  of  p53  and  p38  remained  unchanged 
throughout  the  experiment.  Note  that  contrast  was  increased  for  both 
the  control  and  the  SM-treated  phospho-p38  panels  to  visualize  low 
levels  of  phosphorylated  protein. 


collected  as  described  above.  Protein  (40  /.ig)  from  each  time  point 
was  resolved  by  SDS-PAGE  and  transferred  to  nitrocellulose 
membrane  using  the  iBlot  system  (Invitrogen,  Carlsbad,  CA). 
Western  blotting  was  carried  out  using  standard  washing/incubation 
steps.  Antibodies  (Cell  Signaling  Technology,  Danvers,  MA)  to 
the  following  proteins  were  used,  p38  MAPK  (#9212),  phospho- 
p38  MAPK  (#9211),  p53  (#9282),  and  phospho-p53  (#9284). 


Results 

A  time— course  study  was  performed  in  HaCaT  cells  exposed 
to  SM  to  confirm  activation  of  p53  and  p38,  signaling  proteins 
previously  shown  to  be  activated  in  response  to  SM  (10,  11). 
Western  blot  results  showed  increased  phosphorylation  of  p53 
and  p38  at  20  min  (Figure  1),  suggesting  that  these  proteins 
are  activated  early  after  exposure  to  SM  in  HaCaT  cells. 

SILAC  has  proven  to  be  an  effective  technique  for  identifying 
and  quantifying  hundreds  of  proteins  in  a  single  experiment  and 
is  particularly  suited  for  signaling  studies  in  cell  culture  (5). 
However,  its  feasibility  has  not  been  demonstrated  in  kerati¬ 
nocyte  models  nor  in  studies  involving  chemical  warfare  agents. 
Using  SILAC,  full  incorporation  of  13C15N-labeled  amino  acids 
in  HaCaT  cells  was  accomplished  after  4—5  population  dou¬ 
blings  (Supporting  Information,  Figure  1),  a  standard  time  frame 
for  SILAC  studies.  Titration  of  heavy-labeled  Arg  was  per¬ 
formed  to  determine  the  optimal  concentration  for  minimizing 
light  Arg  and  heavy  proline  (Pro)  production.  We  selected  one- 
third  the  manufacturer’s  recommended  concentration  since  it 
resulted  in  the  lowest  production  levels  of  both  light  Arg  and 
heavy  Pro  (Supporting  Information,  Figure  1),  as  seen  in  our 
earlier  SILAC  studies  (3).  Furthermore,  similar  growth  rates 
and  no  noticeable  differences  in  cellular  morphology  were 
observed  after  three  passages  in  culture  under  dialyzed  FBS 
growth  conditions  (Supporting  Information,  Figure  2),  suggest¬ 
ing  that  SILAC  is  a  suitable  technique  for  use  in  HaCaT  cells. 

Thus,  a  large-scale,  quantitative  analysis  of  the  proteome 
using  SILAC  was  pursued  to  elucidate  additional  cellular 
pathways  and  interactions  activated  in  response  to  SM.  The 
IMAC  method  was  also  employed  due  to  the  importance  of 
protein  phosphorylation  in  many  signaling  cascades  (6).  A 
schematic  of  this  approach  is  shown  in  Figure  2. 

This  method  yielded  9525  total  phosphorylation  sites  on  4760 
phosphopeptides  (Supporting  Information,  Table  1)  with  a  false 
discovery  rate  of  less  than  1%.  A  high  degree  of  correlation 
between  the  experimental  data  set  and  a  predicted  normal 
distribution  of  quantitative  ratios  was  also  observed  (Supporting 
Information,  Figure  3),  and  statistical  analysis  resulted  in  a  mean 
log2-transformed  ratio  of  —0.089  and  a  log?-transformed  SD  of 
0.55.  This  SD  corresponded  to  a  quantitative  ratio  of  1.46, 
indicating  that  phosphopeptide  ratios  in  our  data  set  deviating 
by  at  least  2-fold  fall  outside  two  SDs  from  the  mean. 
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Figure  2.  Flowchart  of  strategy  used  for  large-scale  phosphoproteomic  analysis  of  SM-exposed  cells.  HaCaT  cells  were  cultured  in  the  presence 
of  either  isotopically  “light”  (control,  red)  or  “heavy”  (SM-exposed,  blue)  Arg/Lys.  After  a  20  min  exposure,  lysates  were  collected,  mixed  at  a  1:1 
ratio  by  protein  concentration,  and  resolved  by  one-dimensional  SDS-PAGE.  Gel  fractions  were  excised  and  in-gel-digested  with  trypsin.  After 
IMAC-based  phosphopeptide  enrichment,  samples  were  analyzed  by  high-resolution  LC-MS/MS.  Quantitative  information  was  obtained  by  comparing 
the  intensities  of  peptides  originating  from  the  heavy  and  light  samples. 


After  excluding  redundant  peptides  from  the  raw  data  set 
(Supporting  Information,  Table  1),  our  approach  resulted  in  the 
confident  identification  of  2315  unique  phosphorylation  sites 
(Supporting  Information,  Table  2)  from  a  total  of  790  proteins 
(Supporting  Information,  Table  3).  Peptide  data  were  further 
analyzed  and  summarized  by  phosphorylated  residue  type 
(serine,  threonine,  and  tyrosine)  as  well  as  the  distribution  of 
peptides  according  to  number  of  phosphorylation  sites  (Sup¬ 
porting  Information,  Figure  4). 

Additional  filtering  criteria  were  applied  to  the  list  of  unique 
phosphorylation  sites  to  compile  a  set  of  proteins  containing 
phosphopeptides  altered  by  at  least  2-fold  (increased  or  de¬ 
creased  following  SM  exposure).  In  addition  to  a  2-fold  or 
greater  quantitative  change,  these  criteria  required  that  the 
phosphopeptides  have  a  signal- to-noise  ratio  >5  in  the  MS 
spectrum  (for  both  heavy  and  light  peptides)  and  a  confidence 
score  >85  as  determined  by  the  Vista  algorithm  (18).  Phos¬ 
phopeptides  meeting  these  stringent  criteria  were  combined  and 
sorted  into  proteins,  resulting  in  a  list  of  86  proteins  (Supporting 
Information,  Table  4). 

The  proteins  in  Table  1  of  the  Supporting  Information  were 
further  analyzed  by  Ingenuity  Pathway  Analysis  (IP A,  Ingenuity 
Systems,  Foster  City,  CA).  A  gene  network  was  constructed 


de  novo  from  this  list  based  on  known  molecular  interactions 
from  the  published  scientific  literature,  and  nearly  one-third  of 
these  proteins  (27/86)  was  found  to  associate  through  known 
interactions  (Figure  3). 

Discussion 

The  sequencing  of  the  human  genome  and  the  development 
of  gene  microarray  techniques  have  helped  fuel  a  surge  in 
biotechnology  and  systems  biology  research  efforts  over  the  past 
decade.  Additional  innovations  in  genomic  and  proteomic 
technologies  are  examples  of  how  far  biotechnology  has  come 
in  recent  years.  Consequently,  methods  that  allow  for  studying 
hundreds  or  even  thousands  of  genes  or  proteins  within  an 
organism  are  starting  to  become  routine. 

Likewise,  the  field  of  proteomics  has  been  rapidly  accelerated 
by  the  recent  development  of  fast  and  sensitive  mass  spectrom¬ 
etry  instrumentation.  New  methods  have  also  been  developed 
that  take  advantage  of  stable  isotopes  incorporated  into  the 
proteome,  allowing  for  quantitative  comparisons  between 
experimental  samples.  Consequently,  mass  spectrometry  can 
now  generate  both  qualitative  and  quantitative  data  for  thousands 
of  proteins  in  a  single  experiment  (3). 
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Figure  3.  De  novo  interaction  network  analysis  of  proteins  exhibiting  altered  phosphorylation  levels  in  response  to  SM.  Phosphoproteins  from 
Table  4  of  the  Supporting  Information  were  analyzed  using  IPA  software  to  construct  de  novo  interaction  pathways  as  determined  from  published 
literature.  Phosphoproteins  colored  gray  were  altered  by  at  least  2-fold  in  response  to  SM,  whereas  proteins  colored  white  were  not  detected  in  our 
study  as  exhibiting  any  quantitative  phosphorylation  change.  A  high  confidence  score  of  62  was  determined  for  this  interaction  map  by  IPA.  See 
Table  4  of  the  Supporting  Information  for  IPI  numbers  corresponding  to  each  protein.  Interactions  shown:  acts  on  (closed  arrowhead),  translocates 
to  (open  arrowhead),  inhibits  (T-junction),  inhibits  and  acts  on  (closed  arrowhead  and  T-junction),  and  binding  only  (simple  line).  “Acts  on”  and 
“inhibits”  edges  may  also  include  a  binding  event.  Abbreviations  used  for  interactions  as  determined  by  IPA:  A,  activation;  E,  expression;  I, 
inhibition;  LO,  localization;  M,  modification;  P,  phosphorylation;  PD,  protein— DNA  interaction;  PP,  protein— protein  interaction;  L,  proteolysis; 
RB,  regulation  of  binding;  T,  transcription;  and  TR.  translocation. 


While  proteomic  methods  and  mass  spectrometry  have 
become  powerful  tools  for  studying  various  biological  processes, 
they  have  yet  to  be  fully  applied  within  the  field  of  chemical 
defense.  In  response,  we  have  taken  a  large-scale  quantitative 
proteomic  approach  to  studying  SM-induced  toxicity  using 
stable  isotope  labeling,  phosphopeptide  enrichment,  and  high- 
resolution  mass  spectrometry.  Here,  we  report  a  large-scale 
analysis  of  the  phosphoproteome  in  a  cellular  model  of  SM- 
induced  toxicity. 

This  approach  led  to  the  identification  of  several  proteins 
already  known  to  be  implicated  in  SM-induced  toxicity.  For 
example,  RPS6KA1  (IPI00017305),  a  protein  influencing  NF- 
/cB  signaling  (79),  was  phosphorylated  nearly  3-fold  higher  in 
SM-exposed  cells  in  this  study  (Supporting  Information,  Table 
1).  Earlier  work  by  our  group  also  detected  increased  phospho¬ 
rylation  levels  of  this  protein  following  SM  exposure  in  HaCaT 
cells  (70).  LIG1  (IPI002 1984 1 )  was  previously  reported  to  be 
activated  via  phosphorylation  in  response  to  SM-induced  DNA 
damage  (20).  Consistent  with  these  results,  a  phosphopeptide 
from  LIG1  was  detected  greater  than  2-fold  higher  in  SM-treated 
cells.  A  few  additional  peptides  were  detected  with  reduced 
levels  of  phosphorylation  as  well,  although  the  functional 
significance  of  these  phosphorylation  sites  in  SM-exposed  cells 
remains  to  be  determined. 

While  several  proteins  were  detected  that  were  already  known 
to  be  phosphorylated  upon  SM  exposure,  perhaps  more  interest¬ 
ing  is  the  myriad  of  proteins  detected  that  have  not  been 
previously  characterized  following  SM-induced  toxicity.  Among 
the  proteins  exhibiting  increased  phosphorylation  levels  in  SM- 
treated  cells  were  those  known  to  be  phosphorylated  following 
DNA  damage,  which  occurs  following  SM  exposure  (27).  These 
included  RFC1  (IPI00375358),  SRRM2  (IPI00748920),  and 
SVIL  (IPI000 18370),  each  with  phosphorylation  levels  increas¬ 


ing  4-fold  or  greater.  None  of  these  proteins  have  been 
previously  reported  to  be  implicated  in  SM-induced  toxicity. 

Similarly,  SM  has  been  shown  to  activate  DNA  damage  and 
apoptotic  responses  through  a  variety  of  mechanisms,  including 
the  p53-  and  Bcl-2-mediated  pathways  (7,  22).  TP53BP1 
(IPI00029778)  plays  a  role  in  the  DNA  damage  response  by 
enhancing  p53-mediated  transcriptional  activation  (23),  whereas 
BCLAF1  (IPI00006079)  is  a  Bcl-2-associated  transcription 
factor  involved  in  apoptosis  by  regulating  p53  expression 
following  genotoxic  stress  (24).  Phosphorylation  of  these 
proteins  was  upregulated  by  at  least  5-fold  in  SM-exposed  cells. 

CDK2  (IPI00031681)  can  be  inactivated  when  phosphorylated 
at  Thr-14  or  Tyr-15  (25).  However,  it  has  been  reported  that 
this  Thr-14/Tyr-15  inhibitory  phosphorylation  plays  a  relatively 
minor  role  in  regulating  CDK2  activity  during  normal  cell  cycle 
but  that  this  inhibition  can  also  result  from  DNA  damage,  which 
in  turn  impairs  cell  cycle  transitions  (26).  A  peptide  phospho¬ 
rylated  at  both  residues  was  upregulated  greater  than  3-fold  in 
SM-exposed  cells.  This  suggests  an  inactivation  event  of  CDK2, 
consistent  with  SM-induced  DNA  damage. 

Rho  GTPases  are  known  to  control  a  variety  of  signaling 
pathways  (27).  The  balance  of  GDP  release  and  GTP  hydrolysis 
on  Rho  is  determined  by  guanine  nucleotide  exchange  factors 
(GEFs)  and  GTPase-activating  proteins  (GAPs),  respectively. 
Interestingly,  several  GEFs  and  GAPs  were  identified  in 
this  study.  Phosphorylated  ARHGEF11  (IPI00157442)  and 
ARHGAP21  (IPI00169307)  were  detected  at  2-fold  lower 
concentrations  in  SM-treated  cells,  whereas  phosphorylated 
ARHGAP12  (IPI217418)  concentrations  increased  by  greater 
than  2-fold.  While  the  roles  of  GEFs  and  GAPs  in  SM-induced 
toxicity  are  unclear,  future  experiments  could  identify  the 
downstream  cellular  effects  of  these  proteins  following  SM 
exposure. 
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Changes  in  peptide  phosphorylation  specific  to  SM  exposure 
were  also  characterized  in  this  study.  Indeed,  these  sites  may 
point  to  the  most  interesting  therapeutic  targets  since  the  levels 
of  phosphorylation  are  either  undetected  or  zero  in  unexposed 
cells.  Among  these  are  DPF2  (IPI00023322),  a  predicted 
transcription  factor  involved  in  apoptosis  (28),  and  ZMYND8 
(IPI004 1831 6),  a  protein  known  to  be  phosphorylated  upon 
DNA  damage  (29).  Likewise,  the  phosphorylation  of  these 
proteins  has  not  been  previously  reported  to  be  associated  with 
SM  toxicity. 

An  initial  concern  was  the  absence  in  our  data  set  of  several 
changes  in  phosphoprotein  concentrations  previously  implicated 
in  SM-induced  toxicity,  such  as  p38  and  p53.  This  can  be 
explained  by  various  phosphopeptide  sequences  not  being 
amenable  to  LC-MS/MS  analysis,  low  abundance  and  stoichi¬ 
ometry  of  phosphopeptides,  and  a  limited  sensitivity  of  the 
enrichment  methods.  These  limitations  may  be  overcome  in 
future  studies  via  the  use  of  alternative  proteolytic  enzymes  and 
additional  fractionation  steps.  Nevertheless,  our  results  represent 
a  significant  advancement  in  technological  application  of 
proteomics  to  the  toxicological  mechanisms  following  SM 
exposure.  These  results  also  provide  a  number  of  candidates 
for  follow-up  experiments  and  lead  the  way  for  additional  large- 
scale  proteomics  studies. 

The  added  value  of  identifying  changes  in  specific  phospho¬ 
peptide  concentrations  resulting  from  SM  exposure  allows  for 
the  de  novo  construction  of  SM-specific  protein  interaction 
networks.  Accordingly,  phosphoproteins  found  to  be  altered  by 
at  least  2-fold  in  response  to  SM  and  meeting  additional  filtering 
criteria  (Supporting  Information,  Table  4)  were  analyzed  by  IPA 
software.  The  resulting  network  map  (Figure  3)  shows  a  number 
of  known  interactions  to  include  events  such  as  phosphorylation, 
activation,  inhibition,  expression,  and  protein— protein  and 
protein— DNA  binding.  Interestingly,  Ck2  was  not  found  in  our 
analysis,  while  several  of  its  phosphorylation  substrates  were 
detected.  Additionally,  NF-/rB  and  various  caspases  are  known 
to  respond  to  SM  (30,  31),  and  several  proteins  known  to  be 
phosphorylated  or  otherwise  interact  with  NF-a:B  and  caspases 
were  detected  in  our  study.  These  results  highlight  the  utility 
of  this  method  for  discovering  interaction  partners  for  proteins 
known— as  well  as  unknown— to  be  involved  in  SM-induced 
toxicity.  The  exact  mechanisms  by  which  many  of  these  proteins 
are  phosphorylated  in  response  to  SM  and  their  role  in  toxicity 
will  be  the  topic  of  future  work. 

In  conclusion,  this  study  is  the  first  attempt  to  quantitatively 
characterize  the  phosphoproteome  following  SM  exposure, 
resulting  in  a  list  of  identified  phosphorylation  sites  that  can  be 
used  as  a  starting  point  for  future  validation.  While  some  of 
these  phosphoproteins  could  be  potential  targets  for  therapy, 
these  data  could  also  be  further  mined  using  bioinformatics 
software  to  expand  our  understanding  of  phospho-specific 
interactions  in  SM-exposed  cells.  As  methods  improve  and 
technology  develops,  additional  phosphorylation  studies  can  be 
conducted  to  further  characterize  the  signaling  pathways  altered 
in  response  to  SM  and  other  chemical  warfare  agents. 
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